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A stability problem for the current state of the strong coupling superconductor has been considered 
within the lattice Ginzburg-Landau model. The critical current problem for a thin superconductor 
film is solved within the London limit taking into account the crystal lattice symmetry. The current 
dependence on the order parameter modulus is computed for the superconductor film for various 
coupling parameter magnitudes. The field penetration problem is shown to be described in this case 
by the one- dimensional sine-Gordon equation. The field distribution around the vortex is described 
at the same time by the two-dimensional elliptic sine-Gordon equation. 

I. INTRODUCTION 

An influence of the external magnetic field on the superconducting state near the upper critical field is inten- 
sively studied during the last decade . ft was observed 1 that a strong magnetic field near the upper critical field point 
makes superconducting order parameter fluctuations quasionedimensional and, therefore, strong. Nonrenormalizabil- 
ity of the theory at a strong field leads to impossibility to apply the renormalization group technique. The lattice 
translation symmetry is understood to be crucial for the renormalizability of the theory^. The Ginzburg - Landau 
effective action is shown to acquire the Harper's operator form of the kinetic term that restores a three-dimensional 
nature of fluctuations and renormalizability of the theory. It is important to consider, what physical consequences 
follow from this form of the kinetic term in the Ginzburg-Landau theory. However, a sophisticated mathematical 
structure of this operator makes working with it rather difficult. Here we consider an alternative approach to the 
problem using the quasiclassical approach in the London's limit of extremely small coherence length. 

II. EQUATIONS OF ELECTRODYNAMICS FOR THE NARROW-BAND SUPERCONDUCTOR 

We start from the lattice Ginzburg-Landau functional which was derived ir£ and used for a study of the critical 
behaviour 
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where B =rotA is the magnetic induction, , r = a (T — T c ) /T c , T c is the mean-field approximation superconducting 
transition temperature, e(p) is a 3d-perodic function of the quasimomentum with periods 2ir/ai, 2ir/a,2, 2n/a3. 
Operators V and A(x) commute if we choose the Coulomb-London gauge divA. — 0. This functional can be rewritten 
in the Wannicr 
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representations. The lattice vectors m numerate the lattice cites in the plane perpendicular to the magnetic field. 
Varying Eq. over the vector potential A(x), we obtain the Maxwell equation 
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with the following current density 



f(x)v(-tfiV - — A(x)U(x) + ^(x)v(ifiV - — A(x))V*(x) 
c c 



(5) 



where v(p) =de/dp is the group velocity of the order parameter wave packet. Let us present the complex order 
parameter tp{x) in the form 



?A(x) = i?(x) exp (i</>(x) , (6) 

where R 2 = n s is the density of superconducting electrons. When the coherence length is small in comparison with 
the field penetration length, r is large that is equivalent to a presence of the small parameter in the kinetic term. 
Then the order parameter modulus R can be taken as a spatially uniform solution of Eq. Q : 



t(T)R = -gRl (7) 

almost everywhere. A small parameter before the higher derivative does not guarantee here a complete spatial 
uniformity: this is just a loophole for the vortex core solution to appear. Assuming the order parameter modulus to 
be costant and the phase gradient to be slow we come to the expression for the current in the London limit: 
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where nP s = R 2 ,. Introducing the gauge invariant vector field 
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we can rewrite the Maxwell equation Eq. (J3J) in the form 



c 2e 

— rotrotA = -eRhfi—A). (10) 

47T he 



Eq. 1)11) fl can be easily reduced to the standard London equation 
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with S 2 = mc 2 j (47re 2 i? 2 ) in the continuum approximation limit v(p) = p/m. Eq. I jlOj l will be applied below to 
some classical problems of superconductivity in order to see, what physical consequences follow from an explicit taking 
into account the crystal translation symmetry. 



III. CRITICAL CURRENT IN THE THIN FILM 



We analize the case of the simple tetragonal lattice with the e(p) function of the form 



e(p) = Aj_(2 - cos — cos -j— I + Aiifl - cos — — ). (12) 

h h 11 h 

Let us consider the film of thickness d with the current j parallel to the rr-axis. We assume d << £(T), d << 6(T), 
where £ (T) is the coherence length that secure uniformity across the film thickness respectively of the order parameter 
R and the current density j. Using Eqs. (??, ©, © an d (JHJ, we can write the following free energy minimum conditions 
(neglecting the magnetic field effect as it is usually done in such problems) 
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e x (2e/ (he) A)R + tR + gR 3 = 0. 
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Here 



e x (q x ) = Aj_(l - cos^a), 
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Eliminating v x and e x from Eqs. (|15fl and (?? with a use of the relation 
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e x = A_l 1 - i / I 
we obtain the following relation between j x and R : 
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Introducing dimensionless variables J and / according to formulae 

j x = eaA±Rl J/h, 

R = fRo, 

where Rq = |r| /g is a non-trivial solution of eq. 0, we arrive at the following result 
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where A; = |t| /Aj^. In the limit of hv x / (aA±) —> and, therefore, A — » 0, Eq. I|21|l reduces to the classical form 
for the continuum approximation^: 



J 2 = 2fc/ 4 (l - f 2 ) 
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Two plots for the dimensionless current J dependence on the dimensionless order parameter / are presented in fig. 
1 for the cases of k — > and k = 1. A dependence of the maximum position f m on is depicted in fig. 2. 



IV. FIELD PENETRATION PROBLEM 



Let us return to Eq. (|10|) . It can be used in order to determine the vortex structure outside the core. For the case 
of the kinetic term in the form of Eq. I|12|) we can express Eq. (|1(J[I as following 
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rotrotA = -eR 2 v(eA/ (He)), (23) 



or, in London's gauge divA = 0: 
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Assuming the dispersion law Eq. I|12fl to be valid, we can write Eq. (|24|l for the rr-component: 

V 2 A X - AneRlaA/ (he) sm(2n A x a/^> ) = 0, (25) 
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where $o=2irhc/ (2e) is the Onsager flux quantum. This equation can be reduced for the weak field case to <|llll . 
It is convenient to introduce the variable 

2nA x a/$ Q = X- (26) 

Then Eq. 112511 can be rewritten in the form 

V 2 X -^sin X -0. (27) 

Let us consider the classical one-dimensional field penetration problem basing on Eq. I|25|l . In the one-dimensional 
case it reads 

0-^inX = O. (28) 

Eq. H28fl is the well-known equation for the pendulum (but with different sign) ; the one-dimensional penetration is 
described by the solution^: 

(29) 

where cn{x) is the elliptic cosine; k is the elliptic module. Notice that the field penetration problem for the lattice 
model of a superconductor has a similarity with the fluxon problem in the one-dimensional Josephson contact. 

V. VORTEX SOLUTIONS 

The field distribution around vortices is described by Eq. i|27ll . which is better to rewrite in the form 

9^ + W =p^X- (30) 

This is the elliptic sine-Gordon equation. It is well investigated nowi Screw dislocations in the elasticity theory 
magnetization distribution in easy plane magnets, Berezinskii-Kosterlitz-Thouless transition and other systems can 
be described by this equation. It has a rich variety of solutions including arrays of vortices in the London limit. 
Vortex solutions are determined at the plane, punctured in the vicinity of the vortex cores. The N-vortex solution for 
arbitary positions of vortices can be written in the extreme London limit 5 — > oo as follows 

N 

X=/ n k arctan . (31) 

* v ~ Vk 

Solutions of the nonlinear equation i|3Ufl was found, for instance, inS. A structure of the solitary vortex can be 
analyzed if we note that the radial equation 

<^X_ + 1 dx = J_ sinx (32) 
dr 2 r dr 5 2 

belongs to the third Painleve class^, This equation was thoroughly investigated ini£. Its solutions are shown to have 
logarithmic asymptotics at r — > and oscillate at r — > oo. 

In conclusion, we have considered an effect of the crystal lattice symmetry and boundness of the kinetic term 
operator spectrum on electrodynamic properties of the strong coupling superconductor. 
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FIG. 1: Current - order parameter relation for the thin film. Thick line - the continuum limit; thin line 
case with k = 1. The maximum position gives the critical current value. 
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FIG. 3: Maximum dimensionless current J as a function of k. 



